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d ■ Abstract. The recombination of Frenkel pairs resulting from low energy recoils 

in 3C-SiC has been investigated using first principles and Nudged Elastic Band 
' calculations. Several recombination mechanisms have been obtained, involving direct 

interstitial migration, atoms exchange, or concerted displacements, with activation 
O ! energies ranging from 0.65 eV to 1.84 eV. These results are in agreement with 

experimental activation energies. We have determined the lifetime of the Vsi + Sitc 
Frenkel pair, by computing phonons and the Arrhenius prefactor. The vibrational 
contributions to the free energy barrier have been shown to be negligible in that case. 
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1. Introduction 

Silicon carbide has been extensively studied and used in various applications, due to 
its unique physical, chemical and mechanical properties [Tj. In electronics, SiC is a 
possible replacement for silicon in high temperature, high power and high frequency 
devices. In radioactive environment, possible uses concern fusion reactors, confinement 
matrices or spatial electronics. A good knowledge of the crystalline SiC behaviour 
during and after irradiation is a prerequisite for all these applications. In fact, during 
irradiation, lattice atoms are displaced, resulting in the formation of structural defects 
such as interstitials and vacancies. The damage is accumulating in the material, possibly 
leading to mechanical and electrical properties deterioration, and even amorphization in 
the case of large doses or low temperatures [2]. Also, only ion implantation is a viable 
option to dope SiC-based electronic devices, since dopants have high migration energies 
in silicon carbide, preventing the use of conventional thermal diffusion techniques. The 
accumulation of damage in SiC due to ion implantation has been largely studied [31 H]. 

In order to anneal these defects and recover a good crystalline quality, thermal 
treatments have been applied during or after irradiation. The defects concentration 
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will change according to temperature and irradiation parameters, as a dynamic process 
of defect creation and recombination. To better understand the materials behaviour or 
the crystal recovery, after irradiation or annealing treatement, knowledge of the stability 
and mobility properties of defects is required. In particular, important quantities are 
the activation energies for recombining vacancies and interstitials from the Frenkel pairs 
formed in irradiated materials. Experimentally, it has been shown that distinct recovery 
stages exist as a function of temperature and of the level of damage, for 4H-SiC [5] . In the 
case of 6H-SiC, isochronal and isothermal annealings indicate three recovery stages as a 
function of temperature, with activation energies 0.3±0.15 eV (150-300 K), 1.3±0.25 eV 
(450-550 K), and 1.5 ± 0.3 eV (570-720 K) PEj. This last result is in accordance with 
a study of dynamic annealing in 4H-SiC, yielding an activation energy of 1.3 eV for 
recovery with a temperature ranging from 350 K to 430 K [8j. Finally, in the high 
temperature domain, Itoh et al. have shown that irradiation-induced centers in 3C-SiC 
are fully annealed above 1020 K, with a measured activation energy of 2.2 ± 0.3 eV [9]. 

The available informations are not complete enough to draw a full and 
comprehensible picture of the defects annealing and crystal recovery of silicon carbide. 
It is worth to note that cascade simulations reveal that most of the formed defects due to 
irradiation are Frenkel pairs with very few clustered defects and antisites [10J [11]. Two 
mechanims are then mainly responsible for defect annealings, the first being short range 
recombination of Frenkel pairs, and the second long range migration of point defects. 
The latter has been largely studied, and the activation energies for single point defects 
have been calculated using quantum mechanical approaches [12], [13l [H] . Less knowledge 
has been accumulated for the former process, principally because of the large number 
of possible Frenkel pair configurations. Using classical molecular dynamics, Gao and 
Weber have determined that activation energies for Frenkel pairs recombination range 
from 0.22 eV to 1.6 eV [15]. Another study yields an energy barrier of 1.16 eV for the 
recombination of Vgi and SiSi dumbbell [HH Q2] . These results suffer from the poor 
accuracy of classical potentials for describing transition states. Using more accurate 
Density Functional Theory (DFT) methods and transition states calculation techniques, 
Bockstedte et al. [18] and Rauls [13] have investigated the recombination mechanisms 
and the associated activation energies. However, the few selected configurations have 
been built by associating an interstitial with a vacancy in an arbitrary way. 

Using first principles DFT calculations, we have investigated the recombination of 
Frenkel pairs with short interstitial-vacancy distances in 3C-SiC. The recombination 
mechanism and the associated activation energies have been determined with the 
Nudged Elastic Band (NEB) method [19j [20]. The Frenkel pairs considered in this 
work have been obtained during ab initio molecular dynamics of low energy recoils [21 J . 
These configurations are then supposed to provide a realistic description of Frenkel pairs 
occurring in irradiated SiC, in the case of a low density of defects. After a short report 
of our results on the stability of Frenkel pairs, the calculated mechanisms and energy 
barriers for recombination are described. The lifetime of one Frenkel pair has been 
estimated, taking into account vibrational entropy contributions. Finally we discuss 
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our results in relation with experiments and previous studies. 

2. Methods 

We employed the plane-wave pseudopotential code PWscf included in the Quantum- 
Espresso package [22]. Our calculations have been performed in the framework of 
Density Functional Theory [231 123] and the Generalized Gradient Approximation (GGA- 
PBE) as parametrized by Perdew, Burke and Ernzerof [25]. Vanderbilt ultrasoft 
pseudopotentials have been used for describing electron-ion interactions. We found that 
with these pseudopotentials, a plane wave basis with an energy cutoff of 25 Ry is enough 
to get difference energies converged to 0.01 eV. The silicon carbide lattice parameter is 
a = 4.382 A. Defect configurations are modelled with periodically repeated supercells. 
Single defects and Frenkel pairs stability have been calculated in cells encompassing 216 
and 64 atoms. The Brillouin zone sampling was performed using the T point only for 
large cells, and a 4 x 4 x 4 special k-points mesh for small cells [26] . For recombination 
investigations, smaller cells including 64 or 96 atoms have been considered, as well as 
4x4x4 (64 atoms) or 2 x 2 x 2 (96 atoms) special k-points meshs, in order to keep 
reasonable computational times. In this work, we considered defects in their neutral 
charge state. 

Two different approachs may be used for studying the recombination of Frenkel 
pairs. The first, based on molecular dynamics and Arrhenius plots, allows to determine 
recombination without prior knowledge of the path. Unfortunately, energy barriers for 
Frenkel pairs recombination in silicon carbide are expected to be large. The probability 
to obtain a successful event is then low, and very long molecular dynamics runs at 
high temperature are then required, preventing first principles molecular dynamics 
calculations. It is also not clear whether the harmonic approximation remains valid 
for high temperatures, or whether the right recombination mechanisms are activated. 
The second method, employed in this work, is the determination of the minimum energy 
path for recombination. We used the NEB technique [19] , with three images between 
initial (Frenkel pair) and final (perfect crystal) states. The climbing image algorithm 
was also used to determine precisely the transition state [20]. While such a technique 
allows a precise calculation of the recombination mechanism as well as the associated 
activation energy, a guess of the initial path is required. Non-trivial mechanisms are 
then difficult to obtain. 

3. Frenkel pairs stability 

The Figures [015] shows the considered Frenkel pairs. Structural properties of these 
defects have already been described in a previous publication [27]. In the Tabled] the 
defect formation energies are reported, calculated according to the usual formalism [28J, 
for two different cell sizes. Among all Frenkel pairs, Vc+CSi^Qo] has the lowest formation 
energy. More generally, Frenkel pairs including C vacancies have the lowest formation 
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216 atoms 


64 atoms 




d FP 


Ef 


AE 


Ef 


AE 


V c + C C[ioo] 


0.85 


9.90 


-0.03 


10.69 


-0.17 


V c + CSi[ W0 ] 


0.5 


6.69 


-3.24 


7.06 


-3.80 


Vc + CSi[ 10Q ] 


0.95 


9.96 


0.03 


10.86 


0.00 


*V Si + Si T c 


1.5 


14.08 


-0.44 


14.92 


-1.56 


V Si + Si T c 


0.9 


13.46 


-1.06 


13.26 


-3.22 



Table 1. Calculated formation energy Ef (in eV) of several Frenkel pairs, for two 
cell sizes. cLfp (in o-o) is the interstitial-vacancy distance, and AE is the formation 
energy difference between a Frenkel pair and individual point defects. A larger cell (96 
atoms) has been considered for the configuration * due to the large vacancy-interstitial 
distance. 

energies, as expected regarding the low formation energies of C defects (both vacancy 
and interstitials) compared to Si defects [27] . It would be interesting to compare these 
energies, obtained in the framework of GGA, with previous calculations performed using 
Local Density Approximation (LDA). Unfortunately, the Frenkel pair formation energies 
are often not reported in the available literature. Rauls points to energies of 10-11 eV 
for Vc + CC pairs depending on the distance [T3j . in agreement with our results. 

In the Table [1] is also reported the difference between the formation energy of 
Frenkel pairs and the sum of the formation energies of individual defects. This quantity 
allows to compare the Frenkel pair stability against dissociation. Almost all differences 
are negative, suggesting that Frenkel pairs are more stable than single defects. Also, 
large energy differences are mostly associated with short vacancy-interstitial distances, 
indicating an attractive interaction between vacancies and interstitials [27J. 

One important aspect is the influence of cell size on the stability of Frenkel pairs. 
In the Table [1] is reported the formation energies for cells including 216 or 64 atoms. 
Differences as large as 0.9 eV are obtained. We have observed a similar effect for single 
point defects, with larger formation energies for the carbon vacancy, and interstitials. 
The biggest difference is obtained for silicon interstitials, suggesting that a 64 atoms cell 
is too small to allow the full lattice distortion around the defect. Although present, this 
effect is less pronounced for Frenkel pairs, since there is a better accomodation of strains 
thanks to the vacancy. Using small cells, one may then expect an overestimation of the 
energy barrier for creating the Frenkel pair from the perfect crystal. Nevertheless, it 
is reasonable to assume much smaller differences for the recombination energy barrier, 
since similar lattice distortions are expected for the Frenkel pair configuration and the 
transition state. At most, lattice distortion should be larger for the transition state, 
leading to a slight overestimation of the recombination barriers. 
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Figure 1. Minimum energy path for the recombination of the Frenkel pair Vc+CC[ioo] 
{cLfp = 0.85ao). The insets show the relaxed configurations along the minimum energy 
path. Light (dark) spheres are silicon (carbon) atoms, while the empty circle marks 
the position of the vacancy. 




Reaction path 

Figure 2. Minimum energy path for the recombination of the Frenkel pair Vc + 
C^[ioo] {dpp = 0.5a ). See legend of FigureQ]for further details. 

4. Frenkel pairs recombination 

The Figure [1] shows the recombination path and the energy barrier for the Frenkel pair 
Vc + CC[ioo] (oIfp = 0.85ao). The carbon interstitial migrates along the [100] direction 
until it recombines with the vacancy. The energy barrier is 1.43 eV, which is in the 
upper part of the experimental data. 

The second studied Frenkel pair is Vc + CSi[i 00 ] (dpp = 0.5ao). In that case, 
the recombination path is not direct, and involves the concerted displacement of the 
carbon interstitial together with the silicon atom forming the dumbbell (Figure [2]). 
The calculated energy barrier is 1.24 eV, again in the upper part of the available 
experimental data. This calculated path and energy barrier are in agreement with a 
previous calculation, yielding a similar mechanism and an energy barrier of 1.4 eV for 
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Figure 3. Minimum energy path for the recombination of the Frenkel pair Vc + 
CSi[ W0 ] (cIfp = 0.95a ). See legend of Figured] for further details. 




Reaction path 

Figure 4. Minimum energy path for the recombination of the Frenkel pair Vsi + Sire 
(dpp = 1.5ao). See legend of Figure Q] for further details. 

an equivalent Frenkel pair configuration [IS] , 

The third configuration, Vc + CSi[ W0 ] is close to the previous one, with a larger 
interstitial-vacancy distance {dpp = 0.95a ). Surprinsingly, the enery barrier is much 
lower in this case, with 0.65 eV (Figure [3]). This energy is required for the first step 
of the recombination process, with the transformation of the CSi dumbbell to a CC 
dumbbell, very close to the vacancy. This CC defect then recombines very easily with 
the vacancy. Such a result is in agreement with a previous calculation showing a very 
low recombination energy (0.2 eV) for this last step [18J. 

The Figure H] shows the recombination of the Frenkel pair Vgi + S%tc (d FP = 1.5a ) 
according to an exchange mechanism. The silicon interstitial replaces another silicon 
atom, located between the vacancy and the interstitial original position. This silicon 
atom is filling the vacancy. The whole recombination process is associated with a large 
energy barrier of 1.84 eV. It is likely that another mechanism, more complex and less 
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Figure 5. Minimum energy path for the recombination of the Frenkel pair Vs% + Sitc 
(dpp = 0.9a ). See legend of Figure Q] for further details. 

expensive in energy, occurred for this configuration. In fact, the initial interstitial- 
vacancy distance is large, and one may imagine the formation of a SiSi dumbbell, and 
its subsequent migration toward the vacancy along a longer and easier path. In that case, 
the barrier for the SiSi migration, computed to be 1.4 eV [12], may be the upper energy 
limit for the recombination. We have investigated such a possibility, by computing the 
energy barrier for converting the S%tc shown in the Figure @] to a SiSi[i 00 ] dumbbell. 
The barrier is 1.38 eV, in very close agreement to the calculated energy for the SiSi 
migration [12]. After few steps, the interstitial could recombine with the vacancy along 
the [110] direction, what requires only 0.2 eV [18J. 

The last Frenkel pair configuration, Vsi + Sixc (d FP = 0.9a ), includes the same 
interstitial than the previous one, but with a different local geometry (Figure [5]). 
Here, a trivial recombination mechanism seems the best choice, with a simple 
straight migration of the silicon interstitial through an hexagonal transition state. 
The associated activation energy amounts to 1.03 eV. We also investigated another 
possible recombination mechanism, in which the interstitial and a silicon neighbour are 
exchanged. However, a very large activation energy of 2.37 eV makes this event very 
unlikely. 

5. Average lifetime of a Frenkel pair 

In the framework of harmonic transition state theory, the average lifetime of a Frenkel 
pair is given by r = (1/ Vo)e~ AG//feT '. v$ is the effective attempt frequency, while AG 
is the free energy difference between initial and transition states. Assuming that the 
most important contributions are the defect formation energy and vibrational terms, we 
obtain AG = Ef + AU^ — TASyn,. Using the model of harmonic oscillators, vibrational 
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Frequency (cm ) 

Figure 6. Computed phonons frequencies for the Frenkel pair Vsi+Sirc (bottom) and 
the transition state (top) configurations, both shown in the Figure [5l The inset shows 
the vibrational contributions to the free energy barrier as a function of temperature. 



quantities are easily calculated from the following expressions 
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We have determined these quantities in the case of the recombination of the Frenkel 
pair Vsi + S%tc (d FP = 0.9a ). The phonons frequencies u>i have been calculated in 
the frozen phonon approximation for both the initial and transition states obtained 
in a 64 atoms cell. The figure [6] shows the computed frequencies, as well as the 
vibrational contributions to the free energy barrier as a function of the temperature. 
Our results indicate that vibrational contributions may be neglected for temperatures 
usually considered in experiments. When the temperature is as high as 2000 K, the 
activation energy is increased by 0.1 eV. However, the harmonic approximation is maybe 
not valid in this regime. 

Using the calculated phonon frequencies for the Frenkel pair (FP) and the transition 
state (TS) configurations, the effective attempt frequency is obtained from 



n 



3N 
=1 



CO 



FP 



U 3N-1 TS 



(3) 



The computed value is vq = 298.6 cm' 1 = 8.95 x 10 12 Hz. Using this value and the free 
energy barrier for transition, the average lifetime of the Frenkel pair V$i + S%tc is a f ew 
hours at 300 K, but drops to few fis at 600 K. This result is completely in agreement 
with experiments, showing defects annealing with activation energy around 1.3 eV for 
temperatures higher than 450 K [7]. 
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6. Discussion 

The different recombination mechanisms investigated in this work yield activation 
energies ranging from 0.65 eV to 1.84 eV. The lowest energy is obtained for the annealing 
of a C interstitial, while the largest corresponds to a Si interstitial annealing. This result 
is coherent with the greater mobility of C interstitials in silicon carbide. We found 
that there is a large diversity of possible recombination mechanisms. While simple 
interstitial migration to the vacancy is favoured in some case, lowest recombination 
energies are often associated with more complex mechanisms, involving exchange of 
atoms or concerted displacements. Therefore, simple transition state determination 
methods using constrained paths should be avoided for investigating recombination. 

In this work, we have considered cubic SiC, while most of available experimental 
data concern common hexagonal polytypes such as 4H and 6H. However, the influence 
of the polytypism on the recombination mechanism remains unclear. In a recent study, 
Posselt et al. have compared simple point defects in 4H and 3C structures [29]. They 
have shown that, despite additional stable configurations in hexagonal crystals, structure 
and formation energy are similar for a large majority of defects, due to an equivalent 
local atomic environment. However, one may expect differences in the recombination 
paths, especially for Frenkel pairs with large interstitial- vacancy distances. In fact, 
beyond the first neighbours, hexagonal and cubic structures differ. Nevertheless, it is 
likely that similar recombination mechanisms will occur in both structures, because of 
the equivalent local atomic environment, and that activation energies will be close. 

Experimentally, several temperature stages and activation energies for annealing 
defects have been determined. The lowest measured energy is 0.3 ± 0.15 eV. It may 
be related to the migration of CC dumbbells, requiring an energy of 0.5 eV |12j . and 
subsequent recombination along low energy paths. More probably, it may be related 
to Frenkel pairs recombination with short separation distance. In fact, Bocsktedte 
et al. have calculated the recombination of two Frenkel pairs with energies equal to 
0.2 and 0.4 eV [18]. These specific configurations were not considered in this work 
since we have restricted our calculations to Frenkel pairs obtained from displacement 
energies determinations. Only low-index crystallographic directions are used in such 
calculations, to obtain the extreme values of the displacement energies range. However, 
it is possible that other directions allow the formation of Frenkel pairs with mid-range 
displacement energies but with very low recombination energy barriers. These Frenkel 
pairs will recombine during the first experimental stage. The second and third stages are 
associated with energies of 1.3 ± 0.25 eV and 1.5 ± 0.3 eV, respectively. Our calculated 
barriers fit into these ranges. These stages are then related to annealing of short distance 
Frenkel pairs, or to the long-range migration of interstitials [H]. Finally, an activation 
energy of 2.2 ± 0.3 eV has been reported [9], larger than our computed values. It is 
likely that this stage is associated with the annealing of more complex defects. 
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